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The path-dependent magnetizations (polarizations) and susceptibilities in dipole- and spin-glasses 
are described analytically for the standard protocols of temperature and field variations using the 
phenomenological Landau-type description of multiple metastable states in the nonergodic phases. 
The immediate manifestation of metastable states' multiplicity as staircases of magnetizations (po- 
0^ , larizations) curves in the temperature cycling experiments is explained and described. The obtained 

results are in a reasonable qualitative agreement with the existing experimental data. 
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00 ■ I- INTRODUCTION 

The notorious property of spin- and dipole-glass phases in disordered magnets and ferroelectrics is the dependence 
of a sample state on its previous history, that is on its path in the field-temperature plane to a specific point. Thus 
£h ■ arriving to some point at the plane where such phase exists via different protocols (paths), say, field-cooled (FC) and 
zero-field-cooled (ZFC) ones, we get two different states of a glass system with different net moments, susceptibilities 
etc. And seemingly both these states are quite stable - on the laboratory time scale (hours and days) they show no 
tendency to relax to some unique stable state, see, for example, Refs. Q. In general using the variety of the paths 
to a specific field-temperature point we can get in it a number of such presumably stable states with quite diverse 

S thermodynamic parameters and having different degree of order in the local moment's directions. 
Such behavior strongly implies that the phase space of the glassy system is divided into separate regions by the 
very high free energy barriers so the Arrhenius times for the transitions between these regions are much greater than 
the laboratory times. Now it is not known if this means that the true nonergodicity sets in the spin- and dipole-glass 
phases, i. e. if the free energy barriers diverge in the thermodynamic limit so the system can be confined in the 

1 ^ i ' separate " ergodic component" of phase space for arbitrarily long time in sufficiently large samples. In such a case the 

quasi-static evolution of a glassy system can be described within the notion of " metastable states" - local free energy 
C*~) ■ minima each belonging to the unique ergodic component. 

Such nonergodic picture of the glassy phases with multiple metastable states naturally arises in the mean-field 
spin-glass models Q . In this framework the seemingly nonergodic behavior of random magnets and ferroelectrics can 
be qualitatively reproduced in numerical simulations of the mean- field models [3]-0- It should be noted that this 
^vq agreement with the present day experiments having time spans up to 10 5 sec can fail at longer times as it is quite 
possible that there are actually only finite barriers in real glassy phases and realistic short-range models. Then the 
restoration of the ergodicity at longer times would manifest itself by the slow convergence of the numerous presumably 
metastable states to a true stable one after, say, a month or a year even if the influence of external random fields is 
eliminated in the experiment. Yet till now we have neither definite theoretical nor experimental evidences in favor of 
finite or divergent barriers. 

In spite of this uncertainty in the nature of irreversible phenomena in the glassy phases of disordered magnets and 
ferroelectrics the mean-field paradigm of their nonergodicity can be very useful for (at least qualitative) description of 
these phenomena on the time scales of the order 10 5 — 10 6 sec and, probably, even longer ones. The attractiveness of 
this approach lies in its simplicity and visual clarity. Here the quasi-static evolution of nonergodic system is associated 
with the free energy landscape in which it moves between a number of local minima as this landscape changes under 
field and temperature variations. 

The phenomenological form of such free energy describing the competition between the ferro-phase and the glassy 
one was suggested in Ref. [8|]. This phenomenology is based on the observation [9j that in the mean- field approach 
the mechanism of phase transition in random system can be viewed as a result of the condensation of some collective 
excitations (modes) represented by the delocalized eigenvectors of the matrix of random exchange (in magnets) or 
random matrix of force constants (in random ferroelectrics). In strongly disordered short-range systems the modes 
which can condense first at the glass transition are necessarily those lying at the localization threshold Q, [l(|. 
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According to the random matrix theory [Io| in various short-range random matrix ensembles the modes near this 
threshold have specific sparse fractal structure. It appears that the number of sites participating in them is not 
proportional to the full number of sites N as usual but only to N x , x < 1 . This means that the condensation of 
one such mode results in the appearance of local spontaneous moments in a small fraction of the sample's sites - on 
the sparse fractal set spreading throughout all sample (as the mode is delocalized) and having the fractal dimension 
df = xd < d, d is the sample's dimension. It is argued in Ref. [9| that the condensation of one such mode can not 
stabilize the other fractal modes which do not overlap with the already condensed one. So the condensation of such 
modes will proceed until almost all sites of a sample acquire spontaneous moments. Apparently it needs macroscopic 
number N = N/N x = N 1 ~ d f/ d of modes to cover all sites. As the condensation temperatures of the fractal modes 
are determined by the corresponding eigenvalues which are close to the localization threshold and differ by the order 
1 /N values [ici| the temperature interval in which the sequence of their condensations would take place can be quite 
small - of the order N * l/N = N- d f/ d . 

Thus the phenomenological potential of Ref. 0] is the function of Nq magnetizations (polarizations) m, of these 
fractal modes which describe the prevailing orientations of spins (dipole moments) in different nonoverlapping fractal 
sets of sites consisting of approximately Ni = N x = N df ' d sites. It can have up to 2^° local minima (metastable 
states) which differ by the values and signs of mi. It appears that the transitions between these minima involves 
necessarily the change of sign of at least one m, and this needs to overcome the potential barrier of the order N\. 
Thus we have the divergent in the thermodynamic limit barriers but it can be just the consequence of the mean-field 
character of this phenomenology allowing only the simultaneous upturn of N± spins. 

Meanwhile the less energy consuming paths between the local minima may exist in realistic short-range models. 
Nowadays we have the numerical results for 3d short-range Edwards-Anderson Ising spin-glass certifying that the 
collective excitations in it do have fractal character with df « 2.10 and with the energy cost diminishing with their 
size [ll|, [l2|. The latter point out to the absence of divergent barriers yet one should be cautious extrapolating 
the results of simulations for N = 10 3 to the real systems with N ~ 10 18 . Also one may tentatively suppose that 
sequential upturns of the fractal set of N\ spins on the way between two local minima would create the domain wall 
of a sort between the upturned spins and non- upturned ones with the fractal dimension df — 1 when df > 1. The 
energy cost of such wall is of the order N^ d f~ x ^/ d so we still have the divergent barrier between minima albeit much 
lower than in the mean-field potential. Further numerical simulations on larger samples may shed some light on the 
validity of such arguments for short-range systems. 

With these reservations we may turn to the phenomenological mean-field description of the widely used standard 
protocols (such as FC and ZFC ones) for quasi-static irreversible processes in the spin- and dipole-glass phases. We 
use the basic qualitative feature of the metastable states phenomenology of Ref. [8j for these phases - the existence 
of inclined hysteresis loop filled with magnetization (polarization) curves so the sides of loop are the stability limits 
of metastable states. We show that the difference of the thermodynamic parameters in standard protocols can be 
explained by the trapping of a system in different metastable states resulting from the hysteresis loop's temperature 
evolution. The obtained results allow to explain qualitatively the wealth of experimental data on the temperature 
and field dependencies of net moment in standard protocols [13j - [2l| and in their temperature-cycling modifications 
resulting in the "staircases" of magnetization (polarization) curves [HI], [22| . 



II. METASTABLE STATES AND HYSTERESIS LOOPS 

First we briefly recall the essential results of Ref. Q using the magnetic terminology. According to Ref. § the phe- 
nomenological potential F for the randomly frustrated uniaxial ferromagnet can be expressed via the magnetizations 
nii, i = 1,,Nq, of sparse fractal modes 



F(m) = 




Here \m k ] = Nq ^ m!?, Nq — N^^f ) . The coefficients a > and b > in F are some constants specific for a 

i=l 

given disorder realization, while Tf and r g are linear decreasing functions of temperature T changing their signs at 
temperatures Tf and T g correspondingly also being disorder dependent. 

This F (m) is symmetric under all permutations of rrii and describes the competition of ferromagnetic order with 
the spin glass one represented by the Nq — 1 - dimensional order parameter composed of the independent m, — rrij 
components favoring the antiparallel m t orientations. This potential have up to 2^° metastable states with partially 
ordered to, in cither fully ordered (ferromagnetic) (rrij = mo) or spin-glass phase ([to] = 0). In zero field the ordinary 
thermodynamics predicts the transition between these phases at r/ = r g . We consider here only the spin-glass region 
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T g < 0, Tg < Tf. In a specific crystal Tt is some function of r g which is a linear one at small T g . So we introduce 

t = T g /a, t' = Tf /a = t + ct, t > 0, c < 1. 

The conditions on the parameters to and c imply the presence of only one transition into the glass phase at t = 
in zero field h = after which the crystal always stays in this phase at t < 0. So the homogeneous fully ordered state 
with rrii = mo for all i is only the metastable one at t < 0. Meanwhile in the paraphase at t > it is the only stable 
state. Its magnetization mo obeys the equation 

h/a = t'm + (1 + (3) m%, (3 = b/a. 

and it is stable for 

t + (3 + (3) ml > 

In the inhomogeneous metastable states appearing at t < mi can acquire just two values: m + > and 

m_ = — xm + < 0, (1) 
So they can be characterized by the number N + of m+ values in it or by the parameter 

n = N + /N , < n < 1 

All thermodynamic parameters of a state with a given n can be obtained as functions of n and the parameter x 
defined in Eq. CO). Thus [|| 

m + = V^ti? (x, n) , 
R 0, n) = {l - x + x 2 + 13 [n + (1 - n) x 2 ] }~ 1/2 

and the net magnetization of such state is 

m = M (x,n,t) , (2) 
M (x, n, t) = y—i [n — (1 — n) x] i? (x, n) . 

The expression for the susceptibility 

dm (n, t 7 h) 



X 



dh 



of a state with a given ?i is 



X~ X (a;,n,<)a _1 = t' - t 



l + R{x,n)- 



(1 + x) (2x - 1) (2 - x) + 2/3 [(2a: - 1) n + x 2 (2 - x) (1 - n)] 



(2a; - 1) n + (2 - a;) (1 - n) + 2/3n (1 - n) (1 + x) 
The parameter a: = x{n, t, h) can be found from the equation of state of the form 

H (x, n, t) = h, 

H (x, n, t) I a = (t' -t)M (x, n, t) + y^fix{l - x) R (x, n) 3 
The stability condition for a metastable state in glassy phase reads 

1/2 < x(n,t,h) < 2 



(3) 



(4) 



(5) 



Thus the knowledge of x = x(n, t, h) provides us with the full description of the thermodynamics of metastable states 
and the regions of their existence. Yet to graph the field dependencies of their magnetizations we do not need to 
solve Eq. ((4]) as the couple of Eqs. ([H [4J gives the parametric representations of them with the parameter x varying 
between the stability boundaries of Eq. In the glass phase considered these m = m(n, t, h) are shown in Fig. [1] 
for some values of t, t' and [3. They do not cross one another and fill the interior of inclined hysteresis loop. Such 
loops are often seen in disordered magnets and ferroelectrics, see, for example, Refs. [23j-[25[. 
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As we see the ends of m = m(n,t,h) curves at which metastable states loose their stability may be thought of as 
the upper and lower branches of the loop while the top and the bottom of the loop are bounded by the mo (ft) curve. 
The stability lines x — 2 and x — 0.5 representing these upper and lower branches on (m, ft) plane are described by 
the following parametric equations 

m u (ft, t) : m u = M (2, n, t) , ft = H (2, n, t) , < n < 1 
mi (h, t) :mi = M (0.5, n,t),h= H (0.5, n, t) , < n < 1 

Excluding the parameter n from them we find 

h/a ={t-t)m u --y _ j , 

— h e < h < hf {—m e < m u < m j) ,x = 2: 

h/a = {t-t) mi --y _ j , 

—hf < h < h e (—rrif < m; < m e ) ,£C = 0.5. 

Here the parameters 



h f = H(2.1,t) = aJ — [ i'-ti±4 I . 

7 ^ ' ' ; V 3 + /3 v 3 + /?; 



m/ = M(2,l,t) 



3 + /?' 



/i e = H (0.5, 1, t) = 2a J — ^- ( i' - 4i 1 ~ 



3 + 4/3 V 3 + 4/3/ ' 



m e =M (0.5,1, t) = 24 ' ' 



3 + 4^ 

define the points at which m u and mi join the mo (ft) curve, see Fig. [TJ The useful characteristics of the loop are 
also the field and magnetization of the most disordered state with n = 0.5 at the point in which m(0.5, t, ft) joins the 
lower branch, see Fig. [T] 



ft c = H (0.5, 0.5, t) = a ^^ b (t> tigf ) 
m c = M(0.5,0.5,t) = ^ 15 ^5 ? ". 
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FIG. 2: (color online) The temperature dependence of characteristic fields: hf (solid line), h e (dotted line), hd (dashed line), h c 
(dash-dotted line) for (3 = 0.1, c = 1, to = 3. Directed lines show the temperature and field variations for standard protocols. 



For the temperature region where hf > we can also define the parameter n d (t) of the state having the left stability 
boundary at zero field, see Fig. [TJa). From the equation H (2, n d , t) = we have 



I If 1 1\ 2 2t 



The right boundary field of this state is 

h ln - n(n . A 9(l + /3) (2n d -l) (t'-t)^~t 

h d a = H(0.5,n d ,t) = (7) 

The temperature behavior of these fields is shown in Fig. O Note that for the c < 1 case considered here they all 
stay positive except of hf which can become negative if 

at 

t < tf = (9) 

c - c 

Some basic properties of the glassy phase can be deduced from Fig. [2] Thus above h e the system has unique stable 
state - the homogeneous one. Below hf this state ceases to exist and only some disordered states stay stable here. 



III. STANDARD PROTOCOLS: ZFC AND IR 



The values of the introduced characteristic fields play important role in the determination of the results of the 
standard protocols used in experimental studies of dipole- and spin-glasses. The paths in (h, t) plane for these 
protocols are shown in Fig. [2j Along with the field-cooled (FC) and zero-field-cooled (ZFC) paths there are also 
thcrmo-remanent (TR) and isothermal-remanent (IR) ones which can be seen as the extensions of FC and ZFC via 
the last step - the switching off the field. Mostly the TR and IR protocols are used to obtain the magnetizations - 
the so-called TRM and IRM. 
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Here we should note that many experimental studies do not strictly follow the standard ZFC, IR and TR protocols, 
often after turning on (off) the field the crystal is heated and the parameters recorded during the heating are thought 
to be the ZFC, IR and TR ones [lj. Generally it may not be so, thus obtained data may differ from that of the 
standard procedures and may depend on the temperature at which the cooling ends up. Fortunately, in the glass 
phase such heating procedures give the data identical to the standard ZFC and TR protocols as we shall see later. 
But it is not true for IR protocol. 

Every path (protocol) in Fig. [5] is uniquely defined by some t, h values and now we can find the dependence 
of their resulting thermodynamic parameters on these t, h. Noticing that in real experiments and simulations the 
magnetization follows the outline of the inclined hysteresis loop under quasi-static variations of large-amplitude field 
we may conclude that under the quasi-static conditions: 

a) the system does not leave the metastable state until it becomes unstable; 

b) leaving the unstable state the system goes to the metastable state with the nearest net magnetization. 

These rules seem to be the almost apparent consequences of Langevin dynamics for our potential F (m) and they 
allow to easily reproduce the quantitative features of the above protocols in a glassy phase. It is quite easy to find 
mzFC as cooling at zero field brings crystal to the non-magnetized and most deep [8[ state with n = 0.5. Then 
depending on the magnitude of applied (positive) field the system may stay at this state or enter the lower branch of 
the hysteresis loop or acquire the too value, see Fig. [1] So we have 

!m g , h < h c , 
mi, h c < h < h e , (10) 
m ,h e < h. 

Here m g is the magnetization of n = 0.5 state defined in parametric form as 

m g (h,t) :m g =M (x g , 0.5, t) , h = H (x g , 0.5, t) , (11) 

0.5 < x g < 1. 

Excluding the parameter x g we find the equation for m g in the form 

\h\ < h c (\m g \ < m c ) . 



Also using the above considerations on the evolution of x and n parameters in ZFC process we can find Xzfc from 
Eqs. (HUH}. Thus we have 

X(x g ,0.5,t) , h<h c , 
Xzfc = { x(0-5,n r ,<) , h c <h<h e , (13) 
Xo, h e < h. 



H (0.5, n r , t) = h, 

Xo = a- 1 [f + 3(1+/?) 



m, 



■r 



(14) 



Here n r — n r (h,t) is the parameter of the state to which system arrives at the lower loop's branch. 

The IR process proceeds in the loop interior as shown in Fig. [T]and its description is not more difficult. Evidently 
77iji? (IRM) stays zero at h < h c as then the system does not leave the n = 0.5 state. At higher fields mjR is the 
zero-field magnetization of the state with n = n r which system join at ZFC process, see Eq. (fl4]) . The result of the 
further field growth depends on the hf sign at a given t. If hf > and h > hd then system ends up on the upper 
branch of the loop after switching off the field (Fig. QJ,) while for hf < and h > h e it occurs at the fully ordered 
state (Fig. Hb). Thus 



rrir'd (hd — h) t) (h — h c 
+m u (h = Q) 1? (h - h d 
m r r) (h e — h) t) (h — h c 
+m (h = 0)$(h-h e ),hf < 



_ ■ +m u (h = Q) i? (h - h d ) ,h f >0 , , 

miR ~ \ m r l) (h e - h) t) (h - h c ) 



m r = M (x r ,n r ,t) 
x r (h, t) : H (x r , n r , t) 



to (h = 0) = 



-v 

1+0 ' 
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FIG. 3: (color online) The magnetizations of metastable states m(n,t,h) for = 0.5, c = 0.9, to = !• 



m„ (/i = 0) = M (2, ««,, t) = (3n d - 2) 



3/2 



i'-t 



(16) 



Here is the Heaviside's step function. 

Similarly we get the susceptibility for this protocol 



XlR 



■ X (l,0.5,t)^(/i c -/i) 
+X [x r , n r , t) § (h d - h) i? (h - h c ) 
+Xu(h = 0)ti(h-h d ),h f >0 
X (l,O.5,t)0{h e -h) 
+X (x r , n r , t) d (h e — h) $ (h — h c ) 
I +Xo (h = 0)&(h-h e ),hf<0 

Xu(h = 0)= X (2, n d , t) = 

X (l,0.5,t)=a- 1 (t'-if±f) \ 

Xo (/i = 0) = - (2ai') _1 ■ 



(17) 



(18) 



IV. STANDARD PROTOCOLS: FC AND TR 



To describe the FC and TR parameters we should consider the relative positions of the boundaries of metastable 
magnetizations m = m(n,t, h) defined by Eqs. (J2J SJ) . They are shown for several n in Fig. [3] and we can see that the 
edges of m(n, t, h) sheets form something like an amphitheater (with the infinitesimal steps of order N^ 1 ). So when 
cooling in constant field (in FC protocol) the system comes to the boundary of the fully ordered state the sequence 
of infinitesimal jumps to the states with lower n < 1 begins. This proceeds until the system joins the state with some 
n s which is stable for all t at given h. Then on further cooling the system will be trapped in this state. From Fig. [4] 
showing projections of the x — 2 stability lines of m(n, t, h) sheets on the (h, t) plane we can conclude that this last 
state has the stability line maximum touching the line h = const. So to find the parameter n — n s for this state and 
the value t = t s at which the system becomes trapped in it we should solve the following equations 

n s (h) , t s (h) : H (2, n„ t s ) = h, dH ^ = . (19) 
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-0.4 t -0.2 -4 f -2 



FIG. 4: (color online) The a; = 2 stability lines of metastable states for (a) f3 = 3, c = 0.9, to = 0.1, (c > Co) and (b) /3 = 1, 
c = 0.4, to = 0.7, (c < cq). The states with given n are stable above these lines. 



The solutions to Eqs. fl9]l can be represented as 

2 + z t z (3 + 2/3 - /3z) 



3 ' s 3 2- (l-c)z{3 + 2(3- (3z) 

where z obeys the equation 



(20) 



* 8+ w UJ t( 1 - c ) z ( 3+2 ' 9 -' 3z )- 2 ] = ' ( 21 ) 

< z < 1. 

Once we have this z we can find the magnetization for the n s -state which system joins on cooling at t < t s 

m s (h, t) : m s = M (x, n s , t) , h = H (x, n s , t) , 1 < x < 2. 

Then FC magnetization is 

{to , t s <t,hf < h, 
m u , t s < t < 0, h < hf, (22) 
m s ,t < t s . 

This result is valid when FC path cross the stability lines of metastable states . It is always so if c < cq but for c > Co 
it is only true for fields lower than the maximum reached by hf, cf. Fig. 



^/ max 



2a\ 



t 3 



27 (3 + /3) (c-c )' 



So for c > Co Eq. (j22|) is valid at < h < hf max only. If c > Co, h > /i/ max the system always stays in the fully 
ordered state so 

rriFC = m (23) 
In the same way we get for c < Co or c > Co , h < hf max 

Xo,t s <t,h f < h, 

Xfc = { Xu,t s <t<0,h<h f , (24) 

Xa,t < t s . 
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with 

Xu = X (2, n u , t) , H (2, n u , t) = ft, 2/3 < n u < 1 
Xs = X {x, n s ,t) ,H (x, n s ,t) = ft, 1 < x < 2 

Evidently, for c > cq, ft > ft/ max 

Xfc = Xo (25) 

To get the ititr. (TRM) we should consider the effect of the switching off the field on the FC state. Turning to Fig. 
[T] we conclude that if FC state becomes unstable at ft = then system ends up on the upper branch of the loop, 
otherwise it preserves the state (the n number) down to ft = 0. As 

H(2,n s ,3t s ) = 

the last happens for n s state at t < 3t s while the fully ordered state can stay stable down to ft = at t < tfo when 
c > Co, see Eqs. ([HI HI). Thus we have for c < cq or c > cq, h < hf max 

_ f m u (h = Q),3t s <t <0, , . 

mTR -\m sQ ,t<M s . ^ 

m s0 (h,t) = M (x ,n s ,t) , 
x a (h, t) :H (x ,n a (h) , t) = 0. 



and for c > Co , h > hf , 



_[ffl„(ft = 0),< /o <f<0,(/i / >0), , . 

TR ~ \ m (h = 0),t <t fo ,(h f <0). [Z() 



For xtr we obtain for c < cq or c > cq, ft < ft 



/ max 



and for c > cq , ft > ft / ! 



v _ j X«(ft = 0),3i s <<<0, , 



= i X„(ft = 0),i /o <t<0, (ft/>0), 
X ™ 1 xo(^ = 0) = -(2ai')- 1 ,t<t /0 ,(ft / <0). 1 yj 



The field and temperature dependencies of the obtained thermodynamic parameters for different protocols are shown 
in Figs. ([5l [7]) for some values of c, to and /3. All magnetizations are continuous functions of field and temperature 
but susceptibilities can have jumps. The breaks, jumps and merging points in Figs. (0 [7]) are associated with the 
loop's characteristic fields and the temperatures at which line ft = const, crosses the graphs of these fields in Fig. [21 
i. e. the solutions to equations 

hk(tk) = h, k = c,d,e,f. 

Also we should add to them tfo of Eq.©, t s and 3t s . For c > Co, ft < hf max the relevant tf is the largest root of the 
above equation. There is one more characteristic field ft.3 S for TR process defined as 

3t s (h 3s ) = t. 

Above it ttitr. = m u (h = 0) and does not depend on ft, see Fig[5ji and Eqs. (|261 H6p . 

The jumps of xzfc at h e and t e in Figs. (T5]d, [7]) result from the inhomogeneous approach of m(n, t, ft) to mo when 
n — > 1 and those of xtr, Xir at i = in Fig. [JJare the consequence of — > 2/3 at f — > 0. So the influence of external 
field makes system to be trapped in n = 2/3 state at t — > —0 while at t — » +0 the n = 1 homogeneous state is realized. 
Hence xtr = Xm = Xu{h = 0) = l/at at t -0, sec Eqs. ([TTirTgll^lS^]) . and at t -> +0 x = Xo = l/2at , cf. Eq. 

cm. 

Also mj-R = = rn u (h = 0) vanish at £ — > —0. From Eqs. (JSJ [T5l [551 12"T|) we get at small t 

' -tt 2 1 3/2 
™/«,™«[ 8(3+ ° 2/3) 
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FIG. 5: (color online) Field dependencies of magnetizations (a) and susceptibilities (b) for FC (solid lines), ZFC (dashed lines), 
IR (dotted lines) and TR (dash-dotted lines) protocols for o = 1, j3 = 2, c = 0.9, t = 0.05, t = -0.08. 




FIG. 6: (color online) Temperature dependencies of magnetizations for FC (solid lines), ZFC (dashed lines), IR (dotted lines) 
and TR (dash-dotted lines) protocols for (a)/3 = 3, c = 0.9, to = 0.1, h = 0.02, c > (5/(2 + (3) and (b)/3 = 5, c = 0.5, t = 0.1, 
h = 0.1, c<p/(2 + f3). 




FIG. 7: (color online) Temperature dependencies of susceptibilities for FC (solid lines), ZFC (dashed lines), IR (dotted lines) 
and TR (dash-dotted lines) protocols for (a)/3 = 0.5, c = 0.9, t = 0.1, h = 0.006 and (b)/3 = 2, c = 0.9, to = 0.05, h = 0.006. 

0=1. 
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Here we should note that the experimental observation of the obtained quasi-static magnetic susceptibilities can be 
achieved through the application of a small ac field with sufficiently small frequency. On the temperature and field 
intervals where the system evolves through the succession of states on the upper or lower boundary of the loop it is the 
only way to get the definite susceptibility value as the ordinary quasi-static susceptibility as field derivative of m(h) 
curve does not exist in these regions. Indeed, the increasing of field in the state on the upper boundary brings the 
system inside the loop along the m = m(n, t, h) curve but decreasing of it makes the system follow the loop outline. 
So the left and the right derivatives of the quasi-static m(h) curve are different. Meanwhile the slow and small ac 
field will oscillate along the m = m(n,t,h) curve only and we get the susceptibility defined in the obtained above 
expressions. 

Yet near T g the extremely slow relaxation can prevent the achievement of quasi-static regime at the laboratory 
frequencies. This may explain the absence of experimental data for susceptibilities similar to that of Figs. ([5)3, [7]). 
However the present results on the magnetizations in different protocols conform qualitatively to the experiments. 
The field dependencies of magnetizations (polarizations) in Fig. [5^i agree well with the experimental data of Refs. 

MM- 

The rather unexpected result of the present theory is the prediction that 

lim rriTRjR (h) ^ 

h— *-0 

at t < 0, see Fig. [5Ji. Indeed, from Eqs. (|20l [2T1 l22l [26]) it follows that for h — > the system is trapped in n s = 2/3 
state which is magnetized. Yet the only definite experimental witness in favor of this result we have found in the ttitr 
data of Ref [HI]. Also \fc — Xtr at h — coincide with \ 01 n — 2/3 state, cf. Eqs.(|24l |28|) . According to Eqs.fTTI 
113 [13 Xfc — Xtr at h = merge at the value x(l,0.5,t), see Eq.flTHJ). 

The temperature dependencies of rriFC m Fig. [6] exhibit two essentially different type of behavior, in Fig. [6^, mpc 
grows monotonously at t decreasing while in Fig. [5J3 it has a downward kink at t = tf. This happens due to the 
different c values, c > /3/(2 + /3) in the first case and c < /3/(2 + /3) in the second one. The bot h ty pes oirriFC behavior 
are observed in real spin glasses and ferroelectric relaxors, that of Fig. [5^ was seen in Refs. |l4j . [l6j]-[2(| and that of 
Fig. EJd was registered in Refs. Q, [l5|, [2l[, [24|. In both cases there is a reasonable qualitative agreement with the 
present theoretical results. 

The behavior of m fc an d rnz fc for pure glass phase in Fig. [5] can be compared with that for the case when glass 
phase appears as intermediate one between para- and ferro-phase (see Fig. 12 in Ref. @). The qualitative difference 
is the jump of mzFC i n stronger fields in the last case. It can not appear in pure glass phase so the presence of this 
jump indicates that the glass phase is followed by the ferro-phase at lower temperatures. 

V. TEMPERATURE EVOLUTION AFTER STANDARD PROTOCOLS 

Now we can consider the further evolution of IR and TR states under heating or cooling in zero field. To do 
this we turn to the Fig. [5] where the temperature dependencies of ttitr and to/_r for h = 0.3a are shown along 
with those of several metastable states at h = 0. Here ttitr is the magnetization of the n s ~ 0.98 state below 
t = 3t s (h = 0.3a) = —6.2 and above 3t s it acquires the zero-field m value on the upper branch of the loop (the 
x = 2 stability boundary). Apparently the quasi-static heating will make vtltr to follow its curve but under the 
quasi-static cooling the TR state above 3t s will join some metastable curve as shown in Fig. [5] by the directed lines. 
If after that we heat the system to a higher t and repeat the cooling-heating cycle we find that system evolve in it 
along the metastable state with a lower m(n, h, t). Thus repeating the "cooling - heating for higher t" cycles one can 
get a "staircase" of magnetization's curves, explicitly demonstrating the existence of numerous metastable states in 
spin-glass phase. 

Such experiments with zero-field temperature cycling of TR state were made for iron-nitride fine particles and 
canonical spin-glass, CuMn alloy and for PLZT ceramic relaxor [22j with the results qualitatively similar to that 
of Fig. [51 We may also note that Fig. [5] depicts actually the TRM's for several different fields as for another fields 
they will be composed just of another matestable n s branches joining the x = 2 boundary at another 3t s . The TR 
polarizations for different fields was measured in ferroelectric relaxor CdiNb20-j fl6| and the results are quite similar 
to that of Fig. [8j 

The miR(h = 0.9a) curve in Fig. [SJ crosses at t c (h = 0.9a) = —5.6 < t < td(h = 0.9a) = —3.2 the m(n,h,t) 
lines with 0.5 < n < n,d(td{h = 0.9a)) = 0.82. So in this temperature region the quasi-static zero-field heating or 
cooling of IR state makes system join one of these metastable states (and not follow miR(h = 0.9a) curve). At 
t > td(h = 0.9a) = —3.2 the IR state is on the upper branch of the loop and the situation is similar to that of TR 
state considered above. 

Similarly we can consider the further cooling or heating evolution of ZFC and FC states in the constant field 
achieved in these protocols. Fig. [S^ shows the temperature dependencies of tuzfc and rriFC f° r h = 0.3a along with 
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FIG. 8: (color online) Temperature dependencies of MFcih = 0), Mtr(H = 0.3a), Mm(h = 0.9a) and zero-field magnetizations 
of metastable states for /3 — 5, c = 0.9, to — 1. Directed lines shows the evolution of TR state under cooling-heating cycles in 
zero field. 




FIG. 9: (color online) (a) Temperature dependencies of Mpc(h — 0.3a), MzFc(h — 0.3a) and the magnetizations of several 
metastable states at ft = 0.3a for j3 = 5, c = 0.9, to = 1 (c > (3/(2 + (3)). Directed lines shows the evolution of ZFC state 
under cooling-heating cycles in the field ft — 0.3a. (b) Temperature dependencies of Mfc(H = 0.1a), MzFc(h = 0.1a) and the 
magnetizations of several metastable states at ft = 0.1a for (3 = 5, c = 0.5, to = 0.1 (c < ,0/(2 + (3)). Directed lines show the 
evolution of ZFC and FC states under cooling-heating cycles in the field ft = 0.1a. 

those of several metastable states at the same field. Here mzFC at t < t c (h = 0.3a) = —1.9 is the magnetization of 
the n = 0.5 state in h — 0.3a and at t c (h — 0.3a) < t < t e (h — 0.3a) = —0.46 it is that of the lower loop's branch in 
the same field. Again simple heating makes the system follow the tuzfc curve, while the cooling- heating cycles with 
the subsequently rising upper turning points result in staircase of magnetization curves as directed lines in Fig[5^, 
show. 

The mpcih = 0.3a) at t < t s (h = 0.3a) = —2.1 is the magnetization of the n s — 0.97 state in ft, = 0.3a and at 
t s (h = 0.3a) < t < tf(h — 0.3a) = —1 it is that of the upper loop's branch in the same field. The cooling in the field 
results also in the following the mpc curve but on heating above t s the system joins one of the metastable states 
with n s < n < 1. The last effect could hardly be observed for the case c > (3/(2 + 0) in Fig. [9^. It is much more 
pronounced for c < (3/(2 + (3) as in Fig. ^p. Here the temperature hysteresis cycles can be seen between ttizfc aid 
m FC - 
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We should note that the staircases of metastable magnetization curves starting from TR and ZFC states are also 
obtained in the random magnet having the intermediate ferromagnetic phase between paramagnetic and glassy ones 
[l[ . In this case the boundary of ZFC staircase differs by a steeper rise from that of Fig. [9K in conformity with the 
results of Ref. 0. 

VI. CONCLUSIONS 

The present results allow for the following conclusions: 

1. The reasonable agreement of the present theory with the experiments implies the validity of the suggested in 
Refs. Q mechanism of the phase-space separation in glassy phases. Yet further studies are needed to decide if the 
true nonergodicity sets in them so the description in terms of multiple "metastable states" can be extended beyond 
the present laboratory time scales. 

2. The present theory demonstrates the possible form which the results of future microscopic theory could have to 
describe the irreversible phenomena in random media on (at least) moderate yet macroscopic time scales. 

3. The basic feature of the present approach is the filling of the interior of hysteresis loop with the magnetization 
curves of metastable states. The presence of strongly irreversible behavior only inside such loops seems to be the 
universal property of real dipole- and spin-glasses which should be necessary present in future theoretical developments. 

4. One more essential result is that the states at the loop's sides are on the verge of their stability. So the 
random magnets and ferroelectrics are in the regions of self-organized criticality when they traverse these stability 
boundaries under some field and temperature variations. Underlying physical picture suggests here permanent upturns 
(avalanches) of fractal sets of local spins (dipoles) which can be probably registered in experiments. 

5. It is hard to say if the nonvanishing tofc and ttitr at zero field and the susceptibility jumps are just the 
properties of the present model or the universal features of glassy phases. It seems worthwhile to study these issues 
in real and numerical experiments. 

I gratefully acknowledge useful discussions with V.B. Shirokov, M.P. Ivliev, E.D. Gutlianskii and V.P. Sakhnenko. 
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